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ABSTRACT 

Aims. The helium-peculiar star aCen exhibits interesting line profile variations of elements such as iron, nitrogen and oxygen in 
addition to its well-known extreme helium variability. The objective of this paper is to use new high signal-to-noise, high-resolution 
spectra to perform a quantitative measurement of the helium, iron, nitrogen and oxygen abundances of the star and determine the 
relation of the concentrations of the heavier elements on the surface of the star to the helium concentration and perhaps to the 
magnetic field orientation. 

Methods. Doppler images have been created for the elements helium, iron, nitrogen and oxygen using the programs described in 
earlier papers by Rice and others. An alternative surface abundance mapping code has been used to model the helium line variations 
after our Doppler imaging of certain individual helium lines produced mediocre results. 

Results. Doppler imaging of the helium abundance of a Cen confirms the long-known existence of helium-rich and helium-poor 
hemispheres on the star and we measure a difference of more than two orders of magnitude in helium abundance from one side of 
the star to the other. Helium is overabundant by a factor of about 5 over much of the helium-rich hemisphere. Of particular note is 
our discovery that the helium-poor hemisphere has a very high abundance of 3 He, approximately equal to the 4 He abundance, a Cen 
is therefore a new member of the small group of helium-3 stars and the first well-established magnetic member of the class. For the 
three metals investigated here, there are two strong concentrations of abundance near the equator at longitude roughly 135° consistent 
with the positive magnetic maximum and two somewhat weaker concentrations of abundance near longitude 315° on the equator near 
where the helium concentration is centered and roughly where the negative peak of the magnetic field would be found. Another strong 
concentration is found near the equator at about longitude 45° and this is not explainable in terms of any simple symmetry with the 
helium abundance or the apparent magnetic field main polar locations. 
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j^ • 1. Introduction by iBorra. Landstreet & Thompsonl (Il983l) . Phot ometric spec- 

£r " f - ,. B ro x , n , UD <rtnQ rjTMocoo? troscopic and magnetic field data all vary on the |Norrisj (19681) 

ON ; The magnetic peculiar B (Bp) star a Cen (HR 5378, HD 125823, -^ assumed tQ be the rotation iod of the ^ 

( w> V761 Cen) is arguably one of the most remarkable members ^ ... ,. .. . c ., ,. ... .. c , ,. ., 

• ■ . , _ , b . J , , r: — ; — , „ . — : — n 1 ,,.,-k , Qualitative discussions ol the distribution or helium on the 

r^. . or the Bp class. As noted by Jaschek & Jaschek (1967) and „ e „ , , , vr, — ; — jr~5 r~n J.m4 

017 — 7—; — f— n i,„, i . ,. J . : : \?~~r »t T ^ surface of a Cen have been given by Norris & Baschek (1972) 

Jaschek et al. (1968), low-dispersion spectra show that the MK A \ku 1 ,imo, tu r a ■ 1 u i- 

. . ' — IN , . , r . ., , , , r , and Mihalas (1973). The former proposed a single helium-strong 

C 3 , type of the star, determined primarily by the strengths of the , ., c c „ ,., .f , „ , . , . , 

_V .• • . t^, , , , ^„ t, , , •, 1 , pole on the surface of a Cen while the latter looked at a number 

t-H . He 1 lines, varies between B2V and B8 IV, while the strengths c .. . , .. .. ...... ., ... . , . . 

. . , _ „. ' „. ,. , , ... ' . . b of options including the possibility that the period was twice as 

> of Si 11 and Si in lines observed at higher dispersion remain es- , ,, d • a a f ■ ac *l * > i- w 1- 

, , b r _ _, _ long as the 8.8 period determined from the star s light curve. For 

• y—i , sentially unchanged and correspond to a spectral type of B2. ,, . , odo • jirFT - ! — I JTTvttK 1 j * 

k>( nr; — niT^^n, 1 . . . t-,-^ 1 111. .■ the single wave 8?8 period Mihalas (1973) explored two possi- 

Norns 1968) demonstrated that at the B2 phase the helium lines ,.,... c ., ... ,. ... .. , , .5, . ,, 

rNj 1 ai mi „_ Tr ,,,-, r,oi bihties for the helium distribution: a model with two small caps 

5— i are stronger than in a normal B2 V standard while near B8 the c ... , , ... . ,. .. c ., . c .% a 

rri ' ... , , ,. , r ■ , ,, , of excess helium abundance with an inclination of the star of 45 

. _ w t star displays a pronounced helium deficiency and all except the , .. ... .. ^. C[) , , , , ... 

* J. ,. ,. ,_ /, , ,, • and magnetic obliquity B of 45 , and a second model with a sin- 
strongest He 1 lines disappear. He also found that the He 1 equiv- . c J. v -.u • i- .• c .u r ,ro 
, ■,,,-,, • , „j n ., r^n^A ■ 1 gle cap of excess helium with an inclination of the star of 45 
alent widths, W^, varied on a 8°8 12 + 0?004 period. d fl f 90° 

After the discovery of this extreme helium variability, a Cen £,. ' .. . . „ _ , , 

. ,. . ,. ., £ • t . . j , The magnetic variations of aCen measured by 



was stuaiea oy numerous autnors ror its ngnt ana colour r^ ; tt — TT - ^: 1 i ,r.oo k 1 j 1 

. .... .. . . ,. , .. ° ^ ., . Borra, Landstreet & Thompson (1983) ruled out the two- 

variabihty, line excitation anomalies and line strength varia- ~. — , , , f.,!. . 1 i lr .-- L . „ , 

rr^rin^m c -, . t . . . .• {■ .j . cap or symmetric model of Mihalas (1973) as they found 

tions. Norris (1968) failed to detect a magnetic field stronger ,, ■, , , .• , ~- — <■ T\ . 1 

1 inn ^ 1 . — ^i.-i , £1J i5 the sinusoidal longitudinal magnetic field curve to be nearly 

than 700 G but a relatively weak field was tentatively de- ... ° c .° , A ~ n „ ,» TT 

, , i... ,. D . . — : 1 i , " . 1 , ■ • , A a symmetric with extrema of -470 and +430 G. Maximum He 1 

tected by Wolff & Morrison (1974) and quantitatively confirmed .. . . . 

J 1 j j lne strength occurs near the negative magnetic extremum. 

* Based on observations obtained at the Canada-France-Hawaii ^^V^g the usual Oblique Rotator Model for the star, their 

Telescope (CFHT) which is operated by the National Research Council work then leads to the conclusion that our views of the negative 

of Canada, The Institut National des Sciences de l'Univers of the Centre and positive magnetic poles are very similar so if the inclination 

National de la Recherche Scientifique of France, and the University of is not close to 90°, then (3 must be about 90 °. This then seems 

Hawaii. to be consistent with the ideas presented by iNorris & Baschek 
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(Il972l) and the single cap model of Mihalas (1 19731) where a Cen 
has a single helium-strong cap centered approximately 90° from 
the rotational axis (i.e. near the negative magnetic pole) and the 
star has an inclination of about 45°. 

Considering the star's unique variability and its brightness 
(V = 4.401), it is somewhat surprising that there has been es- 
sentially no detailed investigations of a Cen published since the 
advent of modern solid -state detect o rs in the 1980's until the 
recent paper of Hubrig & Gonzalez (2007). They demonstrate 
that the star's metallic line spectroscopic variability is quite com- 
plex and note that a characteristic of a Cen is that some elements 
such as oxygen have lines where the neutral line strength varies 
in antiphase to the singly ionized lines. They also found weak, 
variable high-excitation emission lines of Si 11, Mn 11 and Fe 11 in 
the spectrum of the star but no emission in the HorHei lines. 
They also noted the split appearence of the Fe 11 and Mn 11 ab- 
sorption line profiles, especially at the helium line strength min- 
imum when the Fe and Mn lines are at maximum strength, and 
comment upon the great difficulty of creating a Doppler image 
given the v sin i of the star is only about 15 km s . 

Our interest in a Cen arose from our investigation of 
the stratification of 3 He and 4 He in several helium-3 stars 
(Bohlender 2005), the appa rent stratification of helium in 
a Cen ( Leone & Lanzafame 119971) . and the suggestion of 
iHunger & Gro ote (1999) that a Cen may be an important tran- 
sitional object where the effects of fractionation of radiatively 
driven winds may become important. We obtained a tentative 
hint of the presence of 3 He in a Cen from a single observa- 
tion of the He 1 /16678 line obtained during the course of the 
helium-3 star investigation and therefore in 2003 we subse- 
quently obtained observations of several He 1 lines throughout 
the 8^8 rotation period of the star. Along with the remarkable 
He 1 line variations it was immediately obvious that the star 
showed pronounced varia tions of metal lines as later reported by 
Hubrig & Gonzalez (2007). Even though the v sin i of a Cen is at 
the extreme low end of what is useful for Doppler imaging and 
our phase resolution is quite coarse, the strong spectrum vari- 
ability encouraged us to try to map the surface abundance dis- 
tribution for the elements He, Fe, N, and O to see if there were 
indications of axial symmetry in these distributions that would 
augment the impressions of the orientation of the magnetic axis 

leane d from the earlier work of lBorra. Landstreet & Thompson! 

1981) . 



2. Spectroscopic data 

Spectra of a Cen were obtained with the Canada-France-Hawaii 
Telescope (CFHT) and the high-resolution Gecko coude spec- 
trograph on nine consecutive nights in 2003 June. Since Gecko's 
mosaic grating is a single-order echellette, the spectra were ob- 
tained at four different grating rotations to obtain 75 to 125A- 
long spectra in regions which included the Hei lines AA4437, 
4471, 6678 and 7281, as well as Ha. The Hei ,17281 line has 
rarely been observed because of strong atmospheric contamina- 
tion in this spectral region but its relative weakness, large isotope 
splitting, and the very high revolving power provided by Gecko 
in this order (R - 158, 000) makes it an ideal diagnostic for 3 He. 
The spectra were processed with IRABJ and included tel- 
luric line removal in the Ha and He 1 ,1728 1 spectral regions us- 




Phase 



Fig. 1. Hipparcos photometry for a Cen plotted using the 
ephemeris given in the text. The photometric minimum in H ap- 
pears to be shifted by 0.05 to 0. 1 cycles relative to this ephemeris 
defined by the u light curve. 



ing the bright, rapidly rotating star a Leo as the telluric stan- 
dard. Typical continuum S/N of the individual reduced spectra is 
about 500 (300 for the ,1728 1 spectral region) and the resolving 
power, as determined by the FWHM of the comparison arc lines 
obtained before and after each grating rotation change, ranges 
from 93, 000 <R< 158,000. 

Phases for each observation were o btained using the most re - 
cent ephemeris determined for a Cen dCatalano & Leona fl996): 



JD(u min ) = 2,442, 808.376 + 8.817744 + 0.000109 x E. 



(1) 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA), Inc., under cooperative agreement 
with the National Science Foundation. 



Th ey improved th e precision of the original period determ ined 
by Norris (1968) by combining data from iNomsl d!971l) and 
iPedersen & Thomsenl d 1977b with new uvby photometry. While 
this ephemeris appears to combine data together to produce a 
well-defined minimum at (p — in m, the helium line strength in- 
dex R shown in their Figure 3 appears to indicate that the helium 
line strength maximum, as well as the vby photometry, precedes 
the m photometric minimum by a small phase increment of per- 
haps 0.05 to 0.10. As shown below, this is significant in that the 
helium maximum for our spectroscopic data appears to precede 
the zero-phase calculated from the above photometric ephemeris 
as well. In Fig. Q]we show that Hipparcos photometry of a Cen 
also shows the same phase shift. This difference in phase appears 
to be real since shifting the Hipparcos photometry to give a min- 
imum at phase zero would require a reduction of the period by 
an amount approximately ten times its quoted uncertainty. There 
has only been about 580 and 1130 stellar rotations between the 
ICatalano & Leo ne (1996) observations, the Hipparcos observa- 
tions and our spectra acquisition so the more recent data should 
only have errors of about 0.007 and 0.014 in phase respectively. 
Based on this discussion we will therefore refer to the ephemeris 
above as identifying the u minimum. 

As noted above, the negative extremum of aCen's magnetic 
field occurs near the phase of maximum helium line strength or 
photometric minimum. 

Table Q] provides a summary of our spectroscopic observa- 
tions. Rotational phases are provided in both decimal and angu- 
lar formats for easy comparison with our models below. A small 
sample of the spectra in a section of a single spectral window are 
presented in Figure [2] 

The number of observations and the uniformity of phase cov- 
erage is much less than ideal for Doppler imaging. However the 
high-resolution of the data and the surprising degree of profile 
shape variations apparent in Figure [2] invites an attempt to ex- 
tract some information about the distribution of elements such as 
He, Fe, N and O over the surface in the hope that the symmetry 
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4412 - 4488A 


R = 


= 93,000 


7238-7321A 


R = 


= 158,000 


HJD Phase 
(2450000+) (Decimal) 


o 


Exp. time 

(s) 


HJD Phase 
(2450000+) (Decimal) 


(°) 


Exp. time 
(s) 



2801.8505 


0.337 


121.3 


600 


2800.7844 


0.216 


77.6 


600 


2801.8579 


0.338 


121.6 


600 


2801.7660 


0.327 


117.6 


600 


2802.8581 


0.451 


162.4 


600 


2802.7766 


0.441 


158.9 


600 


2802.8655 


0.452 


162.7 


600 


2803.7557 


0.553 


198.9 


600 


2803.8453 


0.563 


202.7 


600 


2804.7718 


0.668 


240.4 


600 


2803.8529 


0.564 


203.0 


600 


2805.7523 


0.779 


280.4 


600 


2804.8456 


0.676 


243.5 


600 


2806.8027 


0.898 


323.3 


600 


2804.8531 


0.677 


243.8 


600 


2807.8005 


0.011 


4.0 


600 


2805.8207 


0.787 


283.4 


600 


2808.7530 


0.119 


42.9 


600 


2805.8282 


0.788 


283.7 


600 










2806.8490 


0.904 


325.3 


900 










2807.8477 


0.017 


6.1 


900 










2808.8297 


0.128 


46.2 


600 










2808.8370 


0.129 


46.5 


600 











6506 - 6620A 



R = 95,000 



6592 - 6702A 



R = 96,000 



HJD 

(2450000+) 



Phase 
(Decimal) 



C) 



Exp. time 

(s) 



HJD 

(2450000+) 



Phase 
(Decimal) 



C) 



Exp. time 
(s) 



2801.8222 


0.334 


120.1 


900 


2800.8628 


0.225 


80.9 


900 


2802.8313 


0.448 


161.3 


900 


2801.7982 


0.331 


119.1 


600 


2803.8225 


0.561 


201.8 


900 


2801.8055 


0.332 


119.4 


600 


2804.8245 


0.674 


242.7 


600 


2802.8076 


0.445 


160.4 


600 


2805.7971 


0.785 


282.4 


600 


2802.8148 


0.446 


160.6 


600 


2806.7459 


0.892 


321.1 


900 


2803.7967 


0.558 


200.8 


600 


2807.7535 


0.006 


2.3 


600 


2803.8043 


0.559 


201.1 


600 


2808.8151 


0.127 


45.6 


600 


2804.8019 


0.672 


241.7 


600 










2804.8094 


0.673 


242.1 


600 










2805.7744 


0.782 


281.5 


600 










2805.7820 


0.783 


281.8 


600 










2806.7644 


0.894 


321.9 


600 










2806.7719 


0.895 


322.2 


600 










2807.7725 


0.008 


3.0 


900 










2808.7929 


0.124 


44.7 


600 










2808.8006 


0.125 


45.0 


600 



of distribution will augment what is known about the magnetic 
variability and perhaps indicate the orientation of the magnetic 
axis relative to the rotational axis. 



3. Doppler imaging 

The technique for Doppler imaging of Ap and Bp stars to ob- 
tain a map of the variable abundance of elements over their sur- 
faces has been thoroughly covered in the literature. Examples 
of the c omputing approaches to Doppler imaging are reviewed 
inlRice (1996) and in Pisku nov & Rice! (Il993t) and a compre- 
hensive review of many of the uncertainties in Doppler imaging, 
although focussing in this case mainly on cool stars, is found in 
iRice & Strassm eier (20(J3>. 

A critical issue for Doppler imaging is whether the same 
T e ff can be assumed for the entire surface of the star or whether 
a variable T e g is required in the programming. Fortunately this 
issue has been reviewed in an extensive and careful analysis car- 
ried out by Molnar ( 1974) based on the OAO-2 ultraviolet pho- 
tometry scans. Molnar concluded that neither the photometric 
variability nor the ultraviolet Si line variations were due to a 
variati on in r e ff over the rotation cycle of the star. Krtick a et alj 
(2007) have also demonstrated that the light curve of the helium- 
strong star HD 37776, an object quite comparable to aCen in 
fundamental properties and in having quite pronounced line pro- 
file variations, can be almost completely accounted for by the in- 



homogeneous surface distribution of silicon and helium over the 
star's surface and the bound-free transitions of these elements 
and not surface temperature variations. More recent work on a 
cooler Ap star, HR 7224, also sug gests that surface te mperature 
variations do not exist in Ap stars (Krticka et al. 2009). 

The r e ff and log g values assumed f or aCen by previous 
authors (e.g. Molnar 1 974t iNorrisI [l 97 lb are in the range of 
19,000 K to 20,000 K and logg of about 4.0. We have used 
T e ff = 19,000 and logg = 4.0 for all models discussed in this 
paper. 

A discussion of the adopted rotation axis inclination, i, is 
deferred until the end of this section. 



3.1. Helium 

Ideally, for a helium variable such as a Cen, it would be most 
informative if a Doppler image of the star's helium distribution 
could be obtained with at least moderate resolution. There are 
four He i lines in the four spectroscopic regions given in Table[TJ 
the three singlet lines ^4437.551, 6678.154, and 7281.349, all 
arising from the same 2p 'P lower level, and the triplet A4411 
line and its forbidden component. The latter line is so strong 
and complex as to be unusable for Doppler imaging, given the 
relatively slow rotation of a Cen. As a first very simple step we 
therefore used the 4437A line to obtain an approximate map of 
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irnoqe from He 4437A 




phase = G 



phase=4 



phase=90 



phase= 1 35 



4434 4436 

Wavelength 



Fig. 2. A small sample of the spectra of a Cen discussed in the 
paper. The continuum of each spectrum is plotted at a height cor- 
responding to the phase of the observation calculated using the 
ephemeris given in the text. The relative intensity scale is shown 
by the vertical bar at the lower left. Note the extreme change 
in the strength of the He i A4431 line as well as the pronounced 
profile shape variability of the metal lines. 



the helium distribution. We ignored the complex Stark broaden- 
ing inherent to the line and instead used highly adjusted values 
for radiative, Stark and van der Waals broadening parameters of 
6.0, -3.8 and -6.5 respective ly. Wavelengths and log gf values 
adopted for the helium lines (iKurucz & Belll[l995l) are provided 
in Table [2] 

We show our Doppler image of the distribution of helium 
over the surface of a Cen based on the A4437 line in Fig. [3] 
Where multiple consecutive spectra were obtained we did not 
combine spectra; instead this provided an indication of the noise 
level of the spectra in the plots of line profile fits that follow. 
The fit to the line profiles shown in Fig. [4] is quite good and 
the abundance g eometry is consistent with the single sp ot mod - 
els proposed bv iNorris & Baschekl (I 1 9721) and Mihalas| (1 19731) : 
one hemisphere of the star has an enhanced abundance of he- 
lium while the other hemisphere is very helium deficient. Table|2] 
gives the approximate range of helium abundances determined 
by the Doppler imaging. These are averaged over several adja- 
cent pixels but given the rather poor time sampling should not 
be interpreted too literally. When abundances are averaged over 
the much larger obviously helium-rich and helium-difficient re- 
gions on opposite hemipheres of the star our Doppler imaging 
suggests that the helium abundance likely varies by a factor of 
100 to 150 from one side of the star to the other. 

Interestingly, we discovered some evident problems when 
working with the other helium lines of a Cen to produce a 
Doppler image of the helium surface distribution. These appear 
to go beyond what one might expect to have because of our 
neglect of Stark broadening, non-LTE affects, possible element 
stratification or simply a poor choice of line. In trying to invert 
the profiles of the /16678 line to form an approximate map of 




Fig. 3. An approximate map of the location of the helium con- 
centration on a Cen obtained by a fit to the /14437 helium line. 



the helium distribution it was found that the program required 
an assumption of a stellar radial velocity that was 5 km s _1 more 
positive than was needed for the lines of the other elements as 
discussed in the next section. This seemed to arise because of a 
distortion redward of the center of the line profiles that appears 
roughly between phases 100° and 200°. 

Attempts to reconstruct a helium abundance map using the 
/1728 1 line showed a similar problem as the /16678 line, but even 
more exaggerated. The Doppler image profile fits that should 
match the observed profiles at phases up to about phase 100° 
are confounded by a strong satellite feature about 0.5A to the 
red that persists from about phase 100° to about phase 200° and 
then quickly and fairly completely disappears for the rest of the 
stellar rotation. The attempts by the program to accommodate 
this anomalous feature lead to persistent red-shifted distortions 
of the fitted profiles at other phases (see Fig. |5}. 



profile fit for He 4437A 





Fig. 4. He 1 A4431 line profiles variations of a Cen and the re- 
sults of Doppler imaging. Observed profiles are represented by 
the points while the model fits are solid lines. The phases of the 
observed profiles are marked in degrees on the right side of the 
figure and the wavelength relative to line centre (in A) is indi- 
cated at the bottom. 
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Table 2. Lines used for Doppler imaging of a Cen and the ap- 
proximate range of elemental abundances over the surface of the 
star. 

Wavelength Ion E.P. logg/ Min. Max. 

(A) (eV) " log N/N H log N/N H 



-0.5 C.5 



-0.5 0.5 



4437.551 


He I 


21.242 


-2.034 


-3.9 


-0.2 


6678.154 


He I 


21.242 


0.329 






7281.349 


He I 


21.242 


-0.842 






4432.736 


Nil 


23.415 


0.580 


-4.9 


-2.9 


4447.030 


Nil 


20.409 


0.285 






6610.562 


Nil 


21.600 


0.433 






4414.900 


Oil 


23.442 


0.172 


-4.3 


-1.4 


4416.980 


Oil 


23.419 


-0.077 






4419.596 


FeUI 


8.241 


-2.218 


-5.8 


-2.0 


4431.019 


FeUI 


8.248 


-2.572 








-1 -0.5 0.5 



Fig. 5. As for Fig. g] but for the He 1 Z7281 line profile. 
Conventional Doppler Imaging is not able to reproduce the 
"satellite" feature longward of the line center between phases 
of 100° and 200°. 



There are no known blends from lines of other elements at 
this wavelength but a possible solution to the poor models for 
the /16678 and /17281 lines is presented in Section|4] 



3.2. Oxygen, Nitrogen and Iron 

The lines of ionized oxygen and nitrogen as well as twice ion- 
ized iron listed in Table [2] were found to be the most suitable 
for Doppler imaging from among those available in the two ob- 
served spectroscopic regions. They presented the least problem 
in terms of blending issues, had an adjacent continuum that was 
not significantly distorted by crowding lines or the presence of 
broad wings of helium or hydrogen lines and were sufficiently 
strong to provide good line signal-to-noise. The program used 
for Doppler abundance imaging of Ap stars allows for blends 
in the calculation of the local line profiles but mapping of the 
abundance distribution over the stellar surface is only done for 
one element at a time so blends must either be of a non-variable 
element or blends of lines of the same element (see for example 
the mapping of the 7775A line triplet of oxygen for e UMa in 



Rice, Wehlau, & Holmgren 1997). The dat a for the lin es listed 
in Table |2]come from the VALD database dKupka et al.|[l9 99). 

Doppler images were generated using the lines of Table|2]for 
the elements iron, nitrogen and oxygen. The images were inter- 
compared with one anot her and t hen co mpared with the general 
symmetry discussed by Mihalasl d 19731) and as suggested by the 
helium distribution as generally represented in Fig. [3] 

The images developed from the Fe lines at 4419A and 443 1 A 
are shown as Fig. [6] The observed and model fitting to the line 
profiles for the 4431 A iron line are shown in Fig. [7] The fitting 
to the line profiles for Fe 4419A are very similar. The regions 
of enhanced iron evident on the Doppler images obtained using 
the two Fe lines are well defined and relatively tightly confined 
in area such that the five regions of peak iron abundance pro- 
duce clearly identifiable dips or distortions migrating through 
the line profiles with the phase of the rotation of a Cen. Note 
the very strong similarity of the two images obtained using these 
two lines independently. 

The three nitrogen lines of Table|2]were each used to produce 
separate Doppler images using the same parameters as with the 
iron lines and as shown in Table[3] The line profiles for the three 
nitrogen lines did not show such clear and distinct distortions 
passing through the line profiles with rotational phase as had the 
two iron lines. The Doppler images from each nitrogen line then 
did not show such compact regions of abundance enhancement. 
Fig. [8] shows the average of the Doppler images from the three 
nitrogen lines and it is evident that the image obtained for the 
nitrogen abundance distribution is similar to that for iron in that 
the placement of the abundance peaks is close to the positioning 
of the iron peak abundance patches but as is suggested by the less 
distinct pattern of line profile distortions the abundance patches 
are less extreme and broader in area. 

Similar patterns of abundance distribution were obtained 
from the oxygen lines and, as with the nitrogen distributions, 
the pattern was generally less distinct on the surface of the star 
and the line profiles did not show the clear pattern of migrating 
dips through the line profiles with rotational phase. (Although, as 
with nitrogen, the line profiles exhibited profile distortions that 
changed with phase and that were suggestive of the migrating 
profile distortions seen for the iron lines). The approximate range 
of abundances determined by the Doppler imaging for each ele- 
ment are given in Table [2] 

The geometric parameters used for Doppler imaging a Cen 
are given in Table [3] The choice of 70° for the inclination was 
more convoluted than in other cases of Doppler imaging. The 
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Fig. 6. Doppler images of the distribution of Fe on the surface 
of aCen reconstructed from the profile variations of the 4419A 
and 4431 A lines separately, and an image created from the com- 
bination of the Doppler images of both features. 



usual preliminary search for a minimum in the error of fit of the 
profiles calculated from the best image formed at each of a range 
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Fig. 7. As in Fig.[4lbut for the line profile variations of the iron 
line 4431 A. 



was not highly sensitive to inclination. However, the helium pro- 
files were far better fit with an image formed using an inclination 
of approximately 70°. Exploration of the images produced and 
the profile fits for the Fe, N and O data showed that while the 

average image from N 
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of inclinations to the observed profile data gave the usual broad Fig. 8. A composite Doppler image of the distribution of nitro- 

curve with a minimum error for the iron, nitrogen and oxygen gen as constructed by averaging the maps reconstructed through 

images at smaller values of inclination around 45°. As is also inversion of the line profiles of the nitrogen lines listed in 

usual, the image formed from the line profiles for these elements Table[2] 
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overall fit to the profiles was slightly better at smaller inclina- 
tions, largely due to the better fit to the very bottom of the deeper 
profiles, the inverted bottom of the Fe profiles in particular at the 
two profiles at phase 243° was far better represented in the im- 
ages formed with inclination close to 70°. That distinctive and 
well-defined feature was quite poorly fit at smaller inclinations. 
The choice of microturbulent velocity was made by applying 
the Blackwell method (Blackwell & Shallis 1979) embedded in 
the SPECTRUM program dGrav & Corb allv 1994). The program 
was used with the nitrogen lines and suggested a microturbu- 
lence of 2.5 kms -1 . 



4. Helium Revisited 



Bohlender] d2005l) has used data identical to those presented here 
for a Cen to study the abund ance of 3 He in a nu mber of well- 
known helium-3 stars (Hartoog & Cowlevlll979h . In a Cen the 
peculiar red satellite feature in the /1728 1 line, apparent between 
phases 77° and 240°, is shifted by approximately 0.5A from the 
line center of the latter line. This is intriguingly close to the ob- 
served shift of 0.55A between the 3 He and 4 He components of 
this line (Fred et al. 1951). A similar shift is suggested for the 
/16678 line, consistent with its isotope shift of 0.50A, but this 
feature is less obvious because of the greater line strength. We 
have therefore tried to fit the helium lines by adding a contribu- 
tion from 3 He as described below. 

First, we generated a grid of relatively simple pure hydrogen 
and helium non-LTE model atmospheres with T^ - 19, 000 K, 
logg = 4.0 and a range of 7Y(He/H) from 0.001 to 0.8. We used 
the TLUSTY program (Hubenv & Lanz 1995) and started with 
the published (lLanz & Hubenvll2007l) line-blanketed model at- 
mosphere BG19000g400v2 with r eff = 19, 000 K, logg = 4.0 
and 7V(He/H)= 0. 1 as the input atmosphere for the program. The 
companion spectral line synthesis program SYNSPEC was then 
used to produce grids of specific line and continuum intensities 
for the combined 3 He and 4 He lines for each helium abundance. 
To do so we first had to modify the version of SYNSPEC we used 
(v48 provided by Lanz, priv ate communication) to inco r porate 
broadening parameters of Dimitrijevic & Shal-Brechot (1990) 
for the /17281 line since these were not included in the source 
code. SYNSPEC was also modified so that it would recognize 
the wavelengths of the 3 He components of line profiles and use 
the same broadening parameters as are used for the correspond- 
ing 4 He lines. 

Unfortunately TLUSTY does not permit incorporation of sep- 
arate 3 He and 4 He abundances in the production of a model 
atmosphere. We have therefore assumed that the helium abun- 
dance used for the TLUSTY models is the combined 3 He + 4 He 



Table 3. Adopted parameters for a Cen (HR 5378). 



Parameter 


Adopted Value 


T e /f 


19000 K 


l°gg 


4.0 


vsiru 


15.0 km s _1 


Inclination i 


70° 


Rotation period 


8.817744 days 


Radial velocity 


5.8±0.2kms-' 


Micro turbulence £ 


2.5 km s- 1 


Macro turbulence 


0.0 km s" 1 


Radius 


2.8 solar 



abundance and also that the lower mass of the lighter isotope has 
no affect on the structure of the atmosphere, or on the treatment 
of broadening for the lines. We then included both 3 He and 4 He 
line data in the input line list used by SYNSPEC and artificially 
adjusted the log gf values of the lines to set the relative abun- 
dances of the two isotopes but keep the total helium abundance 
used by SYNSPEC consistent with the abundance used in the 
model atmosphere. For example, in the case of the /17281 line if 
7Y( 3 He + 4 He) = 0.1 in the TLUSTY model and 7Y( 3 He/ 4 He) = 1 
then identical log gf values of (-0.842 - log 2) = -1.143 were 
used for both isotopes in the input line list for SYNSPEC to pro- 
vide /V( 3 He/H) = 0.05 and 7Y( 4 He/H) = 0.05. (The log gf value 
for the 4 He Z7281 line is -0.842.) 

The specific intensities were then used in a disk integration 
program that permits the placement of circular zones of differ- 
ent helium abundances anywhere on the disk to approximately 
model the surface abundance geometry. Each zone is defined by 
an angular radius and its colatitude and longitude relative to a 
point that cr osses the line of sight to the ob server at <f> — (see 
Figure 3 of Bohlender & Landstreetl[l990l for an illustration). 
Bands of different abundances can also be modeled with the ap- 
propriate superposition of two circular patches. At each point of 
the disk integration the local helium abundance is determined 
and specific intensities for the appropriate line and local model 
atmosphere are used to construct the line profile. 

Line profiles were generated for the He i AA4437, 447 1 , 6678 
and 7281 lines. Fitting was conducted on a trial and error basis 
until reasonable fits were found for each line for as simple a con- 
figuration of abundance spots as possible. Our models indicated 
that a large inclination of about 70° was required to reproduce 
the line profiles and it was this result that led us to revisit the in- 
clination derived from Dopper imaging of the metal lines as de- 
scribed in the previous section. Figure [9] shows the resulting fits 
to the /I7281 helium line profiles as a function of phase. While a 
model with just two roughly hemispherical regions with widely 
differing helium abundance produced a reasonably good fit, we 
were able to produce a substantially better fit to the line and the 
red satellite feature at all rotational phases by using three regions 
of different, but uniform, abundances. The illustrated fit is pro- 
duced with a 55° wide helium-rich band with 7V( 4 He/H) = 0.50 
centered at a colatitude of 300° and extending from 25° to 80° 
from what is assumed to be the negative magnetic pole. The 25° 
radius polar region centered at the same colatitude has an ap- 
proximately solar helium abundance of 7V( 4 He/H) = 0.10 The 
rest of the star's surface has 7Y(( 3 He+ 4 He)/H) = 0.004 and 
7V( 3 He/ 4 He) = 1. In other words, a broad helium-rich region 
covers much of one hemisphere of the star, except for small re- 
gion with a normal abundance near the pole, while the other 
hemisphere is extremely helium poor but has a huge relative 
overabundance of 3 He. The 3 He contribution to the line profile 
is readily apparent between phases 0.3 and 0.7. There is not a 
discernible contribution to the helium line profiles from 3 He in 
the helium-rich hemisphere. The broad helium-rich band's co- 
latitude of 300° indicates that it crosses the line of site to the 
observer at approximately the same phase as the Hipparcos pho- 
tometric minimum discussed earlier and is also consistent with 
the phasing of the helium line strength photometry. 

Good fits to the Hei AA4437, 4471 and 6678 lines can be 
produced with very similar models. Minor adjustments to the 
helium abundances are requi red to improve the fits to these lines 
but this is not surprising. Bohlender (2005) has shown that in a 
few of the previously recognized helium-3 stars both isotopes of 
helium show evidence of abundance stratification, with 3 He lo- 
cated higher in the atmosphere than 4 He. Our unstratified mod- 
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Fig. 9. Model fits for the He i /T728 1 line of a Cen in which con- 
tributions from both 3 He and 4 He isotopes are incorporated. The 
continuum of each spectrum is plotted at a height corresponding 
to the phase of the observation calculated using the ephemeris 
given in the text. The relatively intensity scale is shown by the 
vertical bar at the lower left. The helium deficient hemisphere of 
a Cen appears to have a huge overabundance of 3 He. 



els of the A6618 line and especially the A4431 line require some- 
what lower 4 He abundances in the helium-rich regions as well as 
lower 3 He abundances in the helium-weak hemisphere to pro- 
duce satisfactory models. This is exactly what is seen in othe r 
helium-3 stars so we suggest, as did lLeone & Lanzafamd (fl997). 
that both isotopes of helium have a vertical abundance stratifica- 
tion in the atmosphere of a Cen. We will not pursue this question 
in detail until data with better phase resolution are available. 
More significant modifications of TLUSTY and SYNSPEC will 
also be needed to accommodate different vertical stratifications 
of both isotopes of helium in the generation of future (preferably 
fully-blanketed) non-LTE model atmospheres and line profiles. 
Qualitatively, the overall helium abundance geometry agrees 
with our Doppler Imaging result: a Cen is dominated by a 
helium-rich region covering much of one hemisphere of the star. 
The small polar region with somewhat lower abundances is also 
obvious at phase 315° in Fig. [3] However, the most remarkable 
result of the isotope abundance fitting is the fact that the very 
helium deficient hemisphere of the star has a very large relative 
abundance of 3 He. 



5. Discussion 

As was discussed in the introduction to this p aper, the ephemeris 
used here is the photometric ephemeris of ICatalano & Leond 
(1996) where minimum light is at phase but from their Fig. 
3 in their paper where they simultaneously show the light curves 
with the variation of the helium index R it is clear that the he- 
lium maximum preceeds phase by a significant fraction of a 
rotational cycle, perhaps of the order of 0. 10. The crude Doppler 
image of Fig. [3] and the profile fits shown in Fig. [9] are consis- 
tent with this as phase 315° in the image appears to be roughly 
when the maximum helium line strength occurs and the center 
of the helium-rich hemisphere in our simple isotope modeling 
is centered at 300°. The minimum line strength occurs at phases 
of between 120° and 135°. As conjectured bv lMihalasi (1 19731) 
the concentration of helium abundance is at about 90° to the 
rotational pole (i.e. at a latitude of about 0°). The location of 
the maximum of helium abundance appears, from the work of 
iBorra. Landstreet & Thompson! ([1983), to be close to the loca- 
tion of the negative magnetic pole and given the symmetry of 
the magnetic variation, we may assume the positive pole is close 
to latitude 0° near a longitude (phase) of 135°. 

We now want to compare the location of the abundance 
peaks that appear consistently in the multiple images derived 
from the separate lines of iron, oxygen and nitrogen. The loca- 
tion of these peaks for these elements are very similar as noted 
above and shown in Fig. [6] and [8] In these figures we see that 
there are generally five spots. A strong pair evident in the im- 
age at phase 135° and just above the equator. This pair would 
appear to be close to the location of the magnetic positive polar 
region and the helium minimum but significantly these are two of 
the three strongest spots on the surface rather than the weakest. 
Further, they seem to be slightly above the equator rather than on 
it and they straddle the location of the polar region. This bifur- 
cation of the peak abundance into two spots is clearly not some 
artifact of noise in the profile data since as has been pointed out, 
the individual abundance maximum regions for the iron lines are 
so well defined in area that they create clear dips seen migrating 
through the line profiles with phase (Fig. [7j. Matching this pair 
of abundance peaks at phase 315° we see a weaker pair of abun- 
dance peaks near the equator and straddling the location of the 
helium maximum and negative magnetic polar site. 

The short description of the location of four of the five obvi- 
ous abundance concentrations is that there are two stronger con- 
centrations straddling the location of the positive magnetic polar 
region where the helium is weakest and two weaker concentra- 
tions straddling the location of the negative polar region where 
the helium is concentrated. 

This leaves as a puzzle the very evident concentration near 
the equator at phase 45° that produces probably the strongest 
of the features we see migrating through the iron line profiles. 
This feature seems totally unrelated to any of the geometry we 
could associate with a simple dipole-like magnetic field struc- 
ture with poles located as described above or related in any 
way with the apparently simple asymmetric helium distribu- 
tion with its concentration near phase 315°. If we choose to 
assume that abundance concentrations form where the mag- 
netic field lines are perpendicular to the surface of the star, 
then it would appear that a Cen has a more complex field struc- 
ture than just a simple dipole or dipole plus quadrupole pat- 
tern. The somewhat cooler helium-weak star HD 21699 is an- 
other example of a star with a single spot of enhanced he- 
lium abundance on an otherwise extremely helium- weak sur - 
face dBrown. Shore. & Sonnebornlll985t IShavrina et al.ll2010h . 
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It also appears to have a very pec uliar magnetic field geometry 
(iGlagolevskii & Chuntonovl2007l) . 

Inve stigations such as that by 

Bagnulo, Landolfi & Landi Degl'Innocentil ([1999) found 
that while some stars can have their magnetic field modelled 
adequately by a dipole or dipole plus quadrupole field, most of 
their sample required a more complex field. In the particular 
case of the star 53 Cam with its relatively slow ro tation and 
strong magnetic field it was seen dKochukhov et alj ,2004) that 
the field was quite complex. A future attempt at producing a 
Zeeman Doppler image of a Cen will be a challenge, but would 
be of considerable interest since the abundance distribution of 
the lighter elements examined here suggests that the magnetic 
field may deviate considerably from a simple structure. We 
encourage new magnetic field measurements of any kind for the 
the star to greatly improve our understanding of the abundance 
and magnetic field geometry of the star and the relationship 
between the two. 

Once more detailed surface abundance and magnetic field 
geometries are known it would clearly also be of interest to de- 
termine if the relatively large (w 0.06 mag) photometric vari- 
ations observed in the star (Fig. [TJ can be reproduced solely 
from the non-uniform abundances of the star, and in particular 
the extraordinary helium variability. As an aside, we point out 
that despite the very large variation in helium abundance over 
the surface of a Cen there is no discernible variability observed 
in the wings of the Ha profile of the star during a single ro- 
tation. This already suggests that minimal temperature or sur- 
face gravity variations are present over the surface of the star. 
Only in the few A at the core of Ha are profile variations of 
any kind seen. These closely reflect the amplitude and shape of 
the variations seen in the O, N, and Fe lines presented here near 
the same rotation phase and therefore likely arise from small at- 
mospheric structure variations caused by the non-uniform metal 
abundances rather than variations in the helium abundance, tem- 
perature or gravity. 

Our discovery of 3 He in the atmosphere of a Cen adds a new 
member to the small group of helium-3 stars. We believe this 
work also makes it the first firmly established magnetic member 
of this class of objects; aScl a nd HR 7129 are other magnetic 
candidates dHartoog & Cowlevlll979l) that we feel require con- 
firmation of their 3 He content with m odern data. As we have al- 
ready mentioned. iHunger & Groota (119 99) suggested that a Cen 
might represent an important transitional object among the mag- 
netic helium-peculiar stars where the effects of fractionation of 
radiatively driven winds may become important. Our discovery 
that the star is a member of the helium-3 class of stars con- 
firms this idea as does the suspecte d vertical stratification of he- 
lium in its atmosphere reported bv lLeone & Lanzafamd (11997) 
and supported by our data. It is interesting to note the fact that 
a Cen has a r e (j and log g very close to the well-established, but 
non-variable and non-magnetic prototype of this class of stars, 
3 Cen A. Clearly the atmospheres and winds of slowly-rotating 
peculiar stars in this mass range, whether magnetic or not, are 
finely tuned to enable the separation of helium isotope s in their 
atmosphere b y either mass fractionation (IVauclairl 19751) or light- 
induced drift dLeblanc & Michaudll 19931) . 

We have already noted the discovery of weak metallic emis- 
sion lines in a Cen by Hubrig & Gonzalez (2007); the incidence 
of this phenomenon appears to be almost universal among the 
helium-3 stars. 



6. Conclusions 

Doppler images of the abundance distribution of the elements 
helium, iron, oxygen and nitrogen in the helium-variable star 
a Cen appear to be consistent with the general geometry of the 
abundance distribution of helium proposed by iMihalasI (119731) 
where the location of the abundance concentrations is roughly 
at the rotational equator. The helium concentration is centred 
at a lo ngitude between 300° and 315° using the ephemeris of 
Catalano & Leone (1996) and we are able to reproduce abun- 
dance variations most consistently with a relatively large value 
of 70° for the inclination of the star. For the three metallic species 
investigated there are two strong concentrations of abundance 
near the equator at longitude roughly 135°, consistent with the 
positive magnetic maximum, and two somewhat weaker concen- 
trations of abundance near longitude 315° on the equator where 
the helium is concentrated and roughly where the negative peak 
of the magnetic field would be found. Another strong concentra- 
tion is found near the equator at about longitude 45° and this is 
not explainable in terms of any simple symmetry with the helium 
abundance or the apparent magnetic field main polar locations. 
We conclude that the surface magnetic field geometry of a Cen 
is likely more complex that that of a simple dipole or dipole plus 
quadrupole. 

Our detailed analysis of the variations of a number of he- 
lium line profiles show that the helium abundance of a Cen varies 
by a factor of approximately 125 from one hemisphere to the 
other if suspected vertical stratification of helium in the star's 
atmosphere is ignored. A large portion of one hemisphere is he- 
lium rich, except very close to the negative magnetic pole, and 
the other hemisphere is very helium deficient. The latter region 
is also found to have an extreme overabundance of 3 He, with 
7V( 3 He/ 4 He) sj 1. As a result, a Cen is a new member of the 
helium-3 class of stars and the first convincing example of a 
magnetic member of this small group of peculiar objects. 
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